Requirements for contractility in disordered cytoskeletal bundles 
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Actomyosin contractility is essential for biological force generation, and is well understood in 
highly organized structures such as striated muscle. Additionally, actomyosin bundles devoid of 
this organization are known to contract both in vivo and in vitro, which cannot be described by 
standard muscle models. To narrow down the search for possible contraction mechanisms in these 
systems, we investigate their microscopic symmetries. We show that contractile behavior requires 
non-identical motors that generate large enough forces to probe the nonlinear elastic behavior of 
F-actin. This suggests a role for filament buckling in the contraction of these bundles, consistent 
with recent experimental results on reconstituted actomyosin bundles. 



I. INTRODUCTION 

The large-scale motion of living organisms often de- 
pends on their ability to harness the power of nanometer- 
sized molecular motors to generate macroscopic dis- 
placements. For instance, our striated muscles rely on 
clusters — or "thick filaments" — of the molecular motor 
myosin to generate forces. Myosin thick filaments are 
able to slide directionally towards the barbed end of po- 
lar actin filaments (F-actin), and the characteristic or- 
ganization of striated muscles into periodic sarcomeres 
arranged in series allows the transfer of this microscopic 
motion to larger scales [Fig. [lja)] [T]. 

Despite its familiarity, sarcomere-like organization is 
far from a universal feature of contractile actomyosin as- 
semblies. In some instances, partially periodic arrange- 
ments reminiscent of sarcomeres are observed, as in sub- 
cellular contractile bundles known as stress fibers [2]. 
In many other cases, however, no such organization is 
known to exist. Examples include smooth muscle fibers 
[3], transverse arcs @], graded polarity bundles [S], the 
cell cortex [6] and lamellar networks [J]. Sarcomere-like 
contraction is unlikely to apply to these systems, and 
there is no consensus regarding their actual deformation 
mechanism. 

In vitro experiments using purified proteins are use- 
ful for understanding contraction in these systems, and 
have been used to identify the minimum requirements 
of actomyosin contractility since the 1940s [8]. Modern 
attempts using dilute actomyosin gels were not able to 
induce contractility in the presence of actin and myosin 
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alone, although adding the actin cross-linker a-actinin 
did produce observable contraction [9Tfl2]. However, 
a more recent study using denser actomyosin bundles 
shows that F-actin and myosin can induce contractility 
on their own |13j . Unlike in sarcomeres, in these bun- 
dles F-actin lacks polarity ordering and myosins are not 
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Figure 1. Contraction and extension in actomyosin systems, 
(a) Contraction in sarcomeres occurs as motor localized in 
the vicinity of the F-actin pointed ends slide towards F-actin 
barbed ends, (b) A motor located near the F-actin pointed 
ends induces local contraction as in sarcomeres (left), while lo- 
calization near the barbed ends yields extension (right). The 
two effects balance each other in a large class of bundles, 
which we characterize in this paper. 
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aligned in register, and we thus refer to them as "disor- 
dered" . 

Previous theoretical work on disordered actomyosin 
systems include continuum models focused on length 
scales much larger than an individual actin filament (141 — 
[T7] , In these elegant descriptions, contractility is intro- 
duced phenomenologically, which circumvents the ques- 
tion of its emergence from microscopic interactions. Sev- 
eral other studies do however investigate this connection. 
In simulations without sarcomeric organization or cross- 
linkers, thick filaments simply sort F-actin by polarity 
without inducing any overall contraction |18) . To restore 
contractility, several models assume that thick filaments 
tend to dwell at the barbed end of F-actin after slid- 
ing over its whole length [T51 - l2"2"] . or more generally that 
their velocity depends on their position relative to the fil- 
ament [23 . In these models, F-actin tend to have immo- 
bilized motors that transiently act as passive cross-linkers 
at their barbed ends. This essentially introduces a small 
amount of sarcomere-like organization and results in con- 
tractility [18] . However, no direct experimental evidence 
of thick filaments dwelling at the barbed end of F-actin 
is available. 

Here we investigate the possibility of bundle contrac- 
tion in the absence of any sarcomeric organization, in- 
cluding motors dwelling at the filament barbed ends. Af- 
ter introducing a general bundle model in Sec. [TTJ we 
show in Sec. Ill that underlying symmetries between con- 
traction and extension [Fig. |T|(b)] imply that disordered 
bundle contraction requires non-identical motors and a 
non-linear elastic behavior of the filaments. Intuitively, 
non-identical motors induce mechanical frustration in a 
disordered actomyosin bundle, thus generating both con- 
tractile and extensile stresses. The filament non-linear 
behavior then allows the former to deform the bundle 
while resisting the latter, yielding overall contraction. Fi- 
nally, in Sec. |IV| we discuss the limitations of our model 
and propose a minimal model for the contraction of a 
non-sarcomeric bundle. 



II. BUNDLE MODEL 

In order to make general statements about a disor- 
dered bundle of potentially complex internal geometry, 
we develop a detailed description of its mechanics with- 
out resorting to the simplifying mean-field approxima- 
tions widely used in previous studies [TM23] . The main 
assumptions of our model are as follows. First, we as- 
sume that the velocity of a motor only depends on the 
force applied to it by the F-actin to which it is bound. 
Second, the average deformation of a thermally fluctuat- 
ing section of filament in the bundle only depends on the 
force applied longitudinally at its ends. Third, we con- 
sider a stabilized bundle where F-actin polymerization 
and depolymerization do not occur. Finally, our model 
does not include the attachment-detachment dynamics 
of motors and cross-linkers to F-actin. The relevance of 



these choices is discussed in Sec. HVl 

We describe a bundle of arbitrary geometry by subdi- 
viding it into interacting "units" of three types character- 
ized by their lengths, velocities and tensions (Sec. II A I. 
Introducing general notations to describe the physical 
connections between different units, we express the re- 
lationships between their forces and velocities as a func- 



tion of their spatial arrangement (Sec. II B). The central 



physics of bundle mechanics are then described by in- 
troducing the filament force-extension relationships and 
the motor force- velocity relationships (Sec. II C). Con- 
servation equations are then used to derive a compact 
description of the bundle amenable to further discussion 
(Sec. [HP]). 



A. Linkers, filaments and junctions 

Here we consider a single bundle constituted of F-actin, 
thick filaments of myosin molecular motors, and option- 
ally passive actin cross-linkers. F-actin are assumed to 
be aligned in the z direction, and for simplicity we refer 
to the direction of positive (negative) z as the "right" 
("left") in the following. To describe the bundle, we de- 
compose it so as to distinguish three types of "units" 
(Fig. [2}: 

• linker units, representing a whole myosin thick 
filament or an passive actin cross-linker (passive 
cross- linkers are equivalent to immobile motors). 
The total number of linker units in the system is 
denoted as n" . 

• junction units, representing the point of contact 
between a myosin thick filament or passive actin 
cross-linker on the one hand, and an F-actin on the 
other. The total number of junction units in the 
system is denoted as n' . 

• filament units, representing a portion of an F- 
actin comprised between two junction units, or be- 
tween one junction unit and an F-actin free end. 
We do not consider freely floating F-actin (filament 
units with two free ends). The total number of fil- 
ament units in the system is denoted as n. 

We label the filament units by i — 1, . . . , n, and denote 
by fi the tension of filament unit i (fi > for a filament 
unit under extension). We allow filament units to bend 
away from the z-axis while maintaining their overall z 
orientation, and thus introduce filament unit i's contour 
length and end-to-end length as two independent vari- 
ables Li and £j, respectively. Finally, we introduce the 
velocity v\ of the rightmost monomer of filament unit 
i, and the velocity v\ of its leftmost monomer. These 
two velocities are Eulerian velocities; in that respect, one 
might think of junction 1 of Fig. [2] as a bridge and of the 
actin as the water that flows under it. Then the veloci- 
ties v i and t>2 associated with filament units 1 and 2 are 
the velocities of the water just upstream and downstream 
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Figure 2. (Color online) Schematic of a model bundle 
comprising three F-actin, one myosin thick filament bound 
in three different sites (linker 1) and one bound only once 
(linker 2). Colors indicate the division into linker units (red), 
junction units (blue) and filament units (green). Our model 
does not impose any restrictions on the number of junctions 
associated with a linker, whether they involve one or more 
F-actin, or the polarity of the F-actin. 
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Figure 3. (Color online) Summary of notation. A filament 
unit (a), junction unit (b), and linker unit (c) are repre- 
sented, along with the velocities and lengths characterizing 
their state, as well as the forces they are subjected to. 



B. Bundle geometry and topology 

To describe the physical connections between junction 
and filament units, we define the n' x n matrices p and 
A by 



1 if i is the right-hand neighbor of i! 

otherwise 

1 if i is the left-hand neighbor of i' 
otherwise. 



(2a) 
(2b) 



For instance, the bundle represented in Fig.[2]is described 
by 



A = 







Vo 



(3) 

The usefulness of these matrices is illustrated by intro- 
ducing the notation f = (fi, ■ ■ ■ , f n ) for the vector of all 
filament tensions, as well as similar notations L, t, v r , 
v l , v" . The matrix product pf is a vector of length n' 
whose i'th component (pf)i> is the tension of the filament 
unit that is the right-hand neighbor of junction unit i' . 
This neighbor thus exerts a force (pf)i' on junction unit 
i', while its left-hand neighbor exerts — (Xf)i'. The sum 
of these two forces is equal and opposite to the force ap- 
plied to junction unit i' by its linker unit. This last force 
thus reads [(A - p)f\ v [Fig. gb)]. 

Further use of p and A indicates that the velocities 
of the rightmost and leftmost actin monomers involved 
in junction unit i' are (pu')i' an d (Av r )j', respectively 
[Fig. |3jb)]. As junction units are point-like objects, the 
net actin flow in and out of them vanishes and 



from the bridge, as opposed to being the velocities of a 
specific fluid element. The notation is summarized in 
Fig. la). 

We label junction units by i' — 1, . . . , n' , and linker 
units by i" = l,...,n". We introduce the velocity v", 
of linker unit i" , i.e., the velocity of the bridge itself in 
our previous analogy. We choose to work in the refer- 
ence frame where the center-of-mass of all linker units is 
motionless, which reads 

n" 

E <>> = °- (!) 

i" = l 

Note that the velocities v\ and v\ are absolute velocities 
defined in this reference frame, as opposed to being veloc- 
ities relative to the neighboring linker. We delay labeling 
the other forces and velocities involved in junction and 
linker units until the introduction of convenient notations 
in the next section. 



\v r = pv l . (4) 

The reasonings used here can be generalized in the fol- 
lowing way: if X{ IS cL quantity associated with filament 
i, then (px)i> is associated to the right-hand neighbor of 
junction unit i'. For instance, in the two previous para- 
graphs we considered Xi = fi and Xi = v\ or v[, respec- 
tively. Just as this statement relates junction units to 
quantities associated with the neighboring filament units, 
we can conversely relate filament units to quantities as- 
sociated with their junction unit neighbors as follows. If 
x' it is associated with junction unit i' , then (p T x')i is as- 
sociated with the junction unit that is the left neighbor 
of filament unit i if it has one, or is equal to zero if it does 
not (the superscript T denotes the matrix transpose) . A 
similar statement holds for A T . Combining this with the 
fact that each junction unit has exactly one right-hand 
and one left-hand neighbor, we find that (pp T x')ii = x'^ 
always. More generally, 

pp T = \X T = 1, (5) 
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where 1 denotes the identity matrix. We further note 
that filament units may have one or two neighbors, which 
implies 



(p T px) 



if i has a left-hand neighbor 
if it does not. 



(6) 



This means that matrix p T p is a projector onto the sub- 
space of filament units that have a left-hand neighbor. A 
similar statement holds for A T A. This discussion implies 
that n' < n. 

To describe the polarity of the filament units, we in- 
troduce the diagonal n x n matrix EE, where Ha = 1 if 
the pointed end of filament unit i points to the right, and 
IXji = —1 if it points to the left. As an example, Fig. [2] 
has 



n = diag(l, 1,-1, -1,1, 1,1) 



(7) 



A similar diagonal polarity matrix II' is associated with 
junction units, and since neighboring filaments and junc- 
tion units have the same polarity we have 



IT = pUp T = XUX 1 



(8) 



Turning to the linker units, we define the n" x n' matrix 
7 by 



if i' is connected to i 1 ' 
otherwise. 



For instance, Fig. [2] has 



1 = 



1110 
1 



(9) 



(10) 



Therefore, the velocity of the linker unit connected to 
junction unit i' is (-f T v")i' [Fig. [3|b)]. Each junction 
unit is connected to a linker unit, but one linker unit 
can be connected to several junction units, which implies 
n" < n'. 

To obtain the forces associated with the linker units, 
we reason that if i' experiences a force [(A— p)f]i' from its 
linker unit, then the i' exerts an equal and opposite force 
[(p — A)/]j' on the linker unit [Fig. [3jc)] . Force balance 
imposes that the sum of the forces applied on any linker 
unit vanishes, and thus 



E 



[(p-a)/]^ = o, 



connected to : 



or. in vector notation: 



7(p-A)/ = 0. 



(11) 



(12) 



Although this condition comprises n" scalar equations, 
those equations are not all independent. This can be seen 
by considering the mechanical subsystem formed by all 
junction and filament units (but not including the linker 
units). As inertia and friction against the background 



fluid are negligible, the sum of all forces applied to this 
system by junction units must vanish: 



E[(A-p)/], = - j>(p-A)/] 



0. 



(13) 



i"=l 



where the first equality follows from the fact that 7 has 
exactly one element equal to 1 per column and zeros 
everywhere else. As Eq. ( 13 1 is always trivially true, 



Eq. ( 12 ) expresses only n" — 1 linearly independent scalar 



conditions. 



C. Filament elasticity and motor operation 

Having defined notations for all lengths, forces and ve- 
locities in our system (Fig. [3]), as well as enforced velocity 
continuity and force balance conditions, we turn to char- 
acterizing the more substantial physics of the junction 
and filament units. 

Filament units are sections of semiflexible polymers 
shorter than or with a length comparable to their per- 
sistence length. We thus assume that the force required 
to hold a filament unit of given contour length Li in me- 
chanical equilibrium is uniquely determined by specifying 
its end-to-end length li , which defines the force-extension 
relationship F: 



f i = F(l i ,L i 



(14) 



In the case of a thermally fluctuating polymer, li de- 
notes the end-to-end length averaged over thermal fluc- 
tuations. Thermal bending of the filament unit moreover 
implies that is smaller than Lj even when f t — 0. Since 
filament units are shorter than the filament persistence 
length, we expect deformations of this kind ranging from 
zero to rs 20%. In vector notation, we write 



f = F(LL), 



(15) 



Motor operation at junction i' is described by a func- 
tional relationship between the local velocity of the linker 
relative to the F-actin and the force applied to i': 



(pv l 



T 1 
7 V 



),=Vl,{[{\-p)f],}, 



(16) 



where the function V[, a priori depends on the polarity 
11^, of the junction unit. This dependence can be ex- 
plicitly determined by noting that the force-velocity re- 
lationship must not depend on our arbitrary choice of the 
direction of positive z. Reversing this choice is equivalent 
to reversing the sign of all velocities, forces and polari- 
ties. The only way for Eq. ( 16 ) to be invariant under this 
transformation is to write 

(pv l ^v"), = % V V!, {% v [{\ p)f] v ] , (17) 

where function V[, is the force- velocity relationship of mo- 
tor i' , an o priori nonlinear function independent of II'. 
In vector notation, 



pv 



T n 
7 V 



wv [ir (a - P )f] 



(18) 
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D. Dynamical equations for bundle deformation 

To provide a kinematic description of bundle contrac- 
tion and extension, we write the conservation of F-actin 
contour length. The rate of change of the contour length 
of a filament unit is directly related to the velocity at 
which the neighboring junction units slide relative to 
F-actin. If i has a neighboring junction unit on its 
right-hand side, the velocity of the linker unit there is 
(X T j T v")i, and its sliding velocity relative to the actin 
is (X Tr y T v" — v r )i. We rewrite this sliding velocity as 
(X T j T v" — X T Xv r )i, which is equal to the previous ex- 
pression if i has a right-hand neighbor, and to zero if it 
does not. Using a similar reasoning for the left-hand side, 
we obtain 



dL 

lit 



(X T -/ T v" - X T Xv r ) - ( P T 1 T v" - p T pv l ) 



The first (second) term on the right-hand side of this 
equation accounts for actin-linkcr sliding on the right- 
hand (left-hand) side of the filaments units, and vanishes 
for the filaments units that do not have a right-hand (left- 
hand) neighbor junction unit. 

We describe the evolution of the end-to-end length of 
a filament unit in different ways depending on whether 
it has two junction unit neighbors or has a free end. In 
the former case, 



d£i 
dt 



{X T 1 T v") i -{p T ^v") i 



(20) 



where the two terms in the right-hand side are the ab- 
solute velocities of the right and left neighbors of i, re- 
spectively. To describe a filament unit with one free end, 



we first note that it has vanishing tension. If the con- 
tour length of the filament unit is known, its end-to-end 
length is given by its force-extension relation Eq. ( 14 ) 



with fi — 0. Differentiating with respect to time, we 
obtain 



dF dti 

~d£~dt 



dF dLi 
<9L~dT 



= 0. 



Defining the diagonal matrices 





dfi 




dF 


ydLJ 




= § ij 


dL 


/dF" 


\ dfi 




dF 


\9£, 




= § ij 





(21) 

(22a) 
(22b) 



ng\ where Sij is the Kronecker delta, we rewrite Eq. (21 1 as 



d^ 
dt 



dF 

~d£ 



_1 dF dL 
~dL~dt 



(23) 



for a filament unit with a free end [23] ■ Making use of 
Eq. ([g]) and its analog for A T A, we write an equation 
that describes filament units whether they have one or 
two neighboring junction units: 



d£ 
dt 



= (p 1 pX 1 - x 1 Xp 1 )Y v 



-(l-p T P X T X) 



dF 

~dl 



_1 dFdL 
dL dt' 



(24) 



We finally use Eqs . Q , (151 and (181 to eliminate /, 
pv l and Xv r in Eqs. (19) and (24 1. This yields 



dL 

~dt 
di 
d7 



= (p T ~ A T )nV[n'(A - p)F(l, L)] (25a) 
W) 1 M {pTp ~ XTX){pT + AT ) n ' y '[ n '( A - PW, L)] + (p T P X T - X T Xp T h T v>:(i, L), (25b) 



where the nonlinear vector function v"(£, L) is the solution of the linear (in v") system of equations formed by Eq. ([T]) 
and the following vector equation, obtained by combining Eqs. ( 12 ) and ( 15 ), then differentiating with respect to time 
and inserting Eqs. (251 into the result: 



iiP - V%(p t p* t ~ ^VW = 7(p - *) - x T x P T )n'V [n'(A - p )f(i, l) 



(26) 



Equation (26), just like Eq. (12), has only n" — 1 in- 



dependent scalar equations, and supplementing it with 
Eq. |l]) thus results in a complete set of equations for v" . 
Finally, supplementing Eqs. ([!]) , ( p5| ) and ( 26 1 with an 
initial condition [£(t — 0),L(t = Ojjcompletely specifies 
the dynamics of a bundle of arbitrary geometry. 



III. SITUATIONS WITHOUT TELESCOPIC 
DEFORMATION 



We now ask under what conditions contraction occurs. 
We arc interested in bundles much longer than the size 
of any single one of their constituents (F-actin or motor) . 
Significant contraction of such a bundle requires that it 



() 



contracts throughout its length, as opposed to, e.g., con- 
tracting at its extremities while the bulk of the bundle 
retains a constant length. We thus focus on "telescopic 
deformation" , whereby the end-to-end velocity of a bun- 
dle contracting under vanishing external load is propor- 
tional to its length. This constitutes the standard behav- 
ior of contractile actomyosin structures in vivo [TJ I25H27] 
and in vitro [13] . Telescopic deformation is characteris- 
tic of systems formed of a serial arrangement of many 
independently deforming elements, often referred to as 
"contractile units" [151 125) 127] . 

Here we demonstrate two requirements for telescopic 
deformation. In Scc. lIII Al we show that it cannot arise if 
the motors all have identical force- velocity relationships. 
Sec. IIIIBI then tackles situations where motors with dif- 
ferent force- velocity relationships are present. In that 
case, we prove that bundles lacking polarity organiza- 
tion {e.g., sarcomeres) comprised of linearly elastic (i.e., 
rigid) filaments do not undergo telescopic deformation. 
Such rigid filaments units represent situations where the 
filament persistence length is very large {e.g., bundles 
of microtubules and kinesin oligomers), or when the fila- 
ment units themselves are very short {e.g., strongly cross- 
linked bundles, where the spacing between two junctions 
is small). 

To completely determine the bundle dynamics, we need 
to specify an initial condition [£{t — 0),L{t = 0)]. We 
thus choose an arbitrary vector Lq of length n and impose 
L{t = 0) = Lq. To avoid confusion between the effects 
of motor-generated stresses, which are relevant for con- 
tractility, and those of bundle prestress, which are not, 
we consider bundles that are initially stress-free: 



f{t = 0) = F[£(t = 0), L(t = 0)] = 0. 



(27) 



This imposes £{t = 0) = £o, where £q is the solution of 
the equation F{£o, Lq) — 0, and thus completely specifies 
the bundle initial condition. 



and by defining £*{t) as the solution of F[£* (t) , L* (t)] 
i This set manifestly satisfies the initial condi- 
tion [^Q,i ] chosen above. Inserting [£*{t), L*{t)] into 
Eqs. ( [25| , we further verify that these functions satisfies 
the equations of motion, implying that they describe the 
dynamics of the bundle. Using Eqs. ( 18 1 and ( 26 ), we find 
that v" — and pv l = Xv r = H'v* for all times, thus 
confirming that linker units are immobile and that the 
velocities associated with right- and left-pointing actin 
are v* and — v* , respectively. 

In this regime, the maximum relative speed between 
any two actin filament units is 2\v* | irrespective of bundle 
geometry. Defining the contraction velocity of the bun- 
dle as the difference between the velocity of its leftmost 
and rightmost filaments, this implies that the contrac- 
tion velocity cannot exceed the constant 2\v*\. It is thus 
impossible for the contraction velocity to scale linearly 
with bundle length, and telescopic contractility does not 
occur. Instead, filaments are segregated according to po- 
larity, as observed experimentally in Ref . [25] . Note that 
this result does not require the bundle to be disordered. 



Disordered bundles with linearly elastic 
filaments 



In a bundle with arbitrary force- velocity relationships, 
the reasoning of the previous section does not apply, 
and some amount of contraction or extension is gener- 
ally present. We thus ask whether bundles contract on 
average, and find that they do only if filament polari- 
ties are organized across the bundles, or if the filaments 
display nonlinear elastic behavior. 

To prove this statement, we first give a mathematical 
description of bundles devoid of both polarity organiza- 



tion and filament nonlinear elastic behavior (Sec. Ill B 1 1, 



and then use Eqs. (25) and (26) to show that such bun- 



dles do not contract(Sec. IIIB2) 



A. Motors with identical force-velocity 
relationships 

In a bundle where the motors have identical force- 
velocity relationships, the junction units have identical 
spontaneous sliding velocities v* in the absence of ap- 
plied force. Defining v* as the length n' vector with all 
its components equal to v* , we can thus write 

V'{f = 0)=v*. (28) 

We now demonstrate that under this assumption, all 
linker units are immobile throughout the dynamics and 
all right-pointing filaments undergo a uniform translation 
with constant velocity v* , while left-pointing filaments 
translate with — v*. 

We define the set of functions [£* (t) , L* (t)] by 

L* (t) = L + {p T - X T )U'v*t (29) 



1. Mathematical formulation 

Consider an arbitrary bundle, which we denote by 
B. Bundle B is fully characterized by specifying ma- 
trices A, p, II, 7, the force-velocity functions V and 
the initial condition Lq. Therefore, a population of bun- 
dles is fully characterized by specifying the distribution 
V{B) — V{X, p, II, 7, V , L ) of the frequencies at which 
any possible bundle B arises in the population. For a 
population without any polarization organization, V{B) 
must be independent of II, implying in particular that it 
is invariant under polarity reversal: 

r{\, P , n, 7 , V, L ) = v{\ P , -n, 7 , V, L ). (30) 

This relationship clearly does not apply to a population 
of sarcomeres. Indeed, in a sarcomere static cross-linkers 
are restricted to the barbed ends of F-actin while ac- 
tive, mobile motors are found at the pointed ends. Thus 



the polarity-reversed image of a sarcomere is not a sar- 
comere; actually, inverting the polarities of filaments in 
Fig. [lja) results in an extensile, not contractile, struc- 
ture. Assuming in the following that Eq. ( 30 ) holds thus 



excludes sarcomeric contractility from our discussion. 

We further assume that filament units exhibit linear 
elastic behavior, which reads 



dF OF 

-d£ {£ ~ £o)+ dL {L - Lo) > (31a) 



F(£,L 

where dF /dt and dF/dL are constants, i.e. 





fdF\ 




fdF\ 


- 

dt 


\d£ 


= - 
dt 


\dL) 



= 0. 



(31b) 



This assumption is a good description of very stiff fila- 
ments, where t% = Li for any filament unit i. Indeed, 
such filaments can be described by choosing the force- 
extension relationship 



F(£ i ,L i )=K(£ i -L i ), 



(32) 



which satisfies Eqs. (31 1 provided that K is a constant, 



and enforcing the limit K — > +oo. 



2. Proof of the property 

To determine whether a bundle B — 
{A,p,n, 7, V',Lq} is contractile or extensile, we 
imagine labeling its leftmost and rightmost points, and 
ask whether the distance C B between these two labels 
tends to increase or decrease with time. To calculate 
C B , we choose a path along the bundle's linker and 
filament units goingfrom the left label to the right label 



as pictured in Fig. |4| We define ef as equal to 1 if the 
path considered crosses filament unit i from left to right, 
to —1 if it crosses it from right to left, and to otherwise 
(see the caption of Fig. [4]). The contraction velocity of 
the bundle can then be defined as: 



dC b 
dt 



df B 

^ 1 dt ■ 



(33) 



We use the following notation to refer to the solution of 
the equations of motion for bundle B: 



L it) = Lo + AL (t) 
£ B {t)=£ + A£ B {t). 



(34a) 
(34b) 



We now introduce bundle B as the polarity-reversed 
image of B, i.e., B — {A, p, —II, 7, V, Lp}. Substitut- 



ing Eqs. (31 1 into Eqs. (25) and (26), we find that the 
dynamics of B satisfies 

AL B (t) = -AL B (t) 
A£ B {t) = -A£ B {t). 



(35a) 
(35b) 
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Figure 4. (Color online) Example of a path (red dotted line), 
as discussed in Sec. |III B| The leftmost and rightmost points 
of the bundle are labelled by gray circles. In this example 
ef — ef = ef 4 = ef Q = — ef = 1 a nd all other ef s are 
equal to zero. As a consequence, Eq. (33 1 reads d£ s /dt = 
d(ii +if -If + £f 4 +^f )/dt 



Combining this with Eqs. (33) and (34) while using the 
path ef = ef to assess the contraction of B, we find 



dC b 
~dT 



dC b 
dt 



(36) 



We finally calculate the average contraction velocity 
over a population of bundles as 



dC 
~d~t 



Em 



dC B 
dt 



(37) 



where the sum runs over all possible bundles. Reorganiz- 
ing this sum, we find 



dC 
~dt 



IE 



V(B) 



dC B 
dt 



dC 1 



- dC B 



dt 



dC B 
dt 



= 0, (38) 



where Eqs. (30) and (36) are used to derive the second 



and third equalities, respectively. Eq. (38) demonstrates 
that bundles without polarity organization or nonlin- 
ear elastic behavior do not contract or extend on av- 
erage. This result does not depend on bundle struc- 
ture or the form of the motor force-velocity relation- 
ships, and can easily be generalized to bundles pinned 
to a rigid substrate, or to include friction of the linker 
or filament units with the solvent. Mean-field model- 
ing of dilute actomyosin gels with rigid filaments suggest 
that this symmetry-based reasoning could have a three- 
dimensional counterpart [23] . However, geometrical non- 
linearities in two or more dimensions can take on the role 
played by elastic nonlinearities in one-dimensional bun- 
dles, thus enabling contraction in disordered networks of 
rigid filaments |29 . 



IV. DISCUSSION 

In this paper we consider bundles of sliding motors 
and filaments with arbitrary geometries and motor force- 



velocity relationships and show that their contractility 
requires 

1. non-identical motors 

2. and 

(a) either polarity organization 

(b) or nonlinear elastic response of the filaments. 

While our model is framed in term of F-actin and myosin 
for clarity, our results are much more general and could 
equally apply to bundles comprised of other types of mo- 
tors and filaments (e.g., kinesins and microtubules). Our 
description includes as a special case the well-understood 
contractility of striated muscle sarcomeres [Fig. [ija)]. 
Their architecture includes both identical myosin thick 
filaments and passive cross-linkers (which are mathemat- 
ically equivalent to motors with velocity zero), thus sat- 
isfying condition (1). They moreover have a distinctive 
polarity organization, and thereby fulfill condition (2a). 
Similarly, bundles with motors whose velocities depend 
on their position relative to the filaments [T9H23] gener- 
ically break polarity-reversal symmetry, which results in 
polarity organization. 

Besides establishing the requirements for contractility, 
the formalism presented here can describe the dynam- 
ics of a wide range of contractile bundles. For instance, 
straightforward numerical simulations of Eqs. (25) and 



(a) 



EE 



(b) 



(26) could be used to describe bundle deformation as a 



function of the initial arrangement of the filament and 
motors. Comparing these predictions to experimental 
observations while varying these initial parameters could 
yield insight into the architecture of the bundles, which 
is currently lacking. Although such a study is beyond the 
scope of our work, a simplified version of our formalism 
still successfully predicts the onset of their contraction 

While the model used here is designed to describe a 
large class of contractile bundles, some of our assump- 
tions are especially appropriate to describe the reconsti- 
tuted bundles of Ref. [13] . F-actin is phalloidin-stabilized 
in this system, implying that no actin polymerization- 
depolymerization takes place; myosin thick filaments are 
much shorter (~ 300 nm) and thicker (~ 50 nm) than 
F-actin (~ 5/xm and ~ 5nm, respectively), justifying 
our assumption that they behave as rigid objects; and 
myosin thick filaments do not detach from the bundle on 
the time scales relevant for contraction, suggesting that 
the filament-motor attachment-detachment dynamics are 
inessential to contractility. Indeed, we show elsewhere 
that attachment-detachment can limit contractility un- 
der low myosin conditions (30) . 

Going beyond the specifics of the system studied in 
Ref. [15] . it is interesting to discuss bundles where 
this attachment-detachment dynamics is not negligible. 
In the simplest such situation, motors undergo attach- 
ment and detachment at a constant rate, which tends 
to randomize their distribution in a filament polarity- 
independent manner. This fails to break the polarity- 



strong weak 
contraction expansion 



weaker stronger 
contraction expansion 



Figure 5. The Fenn effect impedes bundle contractility, 
(a) Simple contractile bundle comprised of fast (light) and 
slow (dark) motors. Double-headed arrows indicate local con- 
traction or expansion due to the movement of the motors 
relative to the filaments [Fig. [ljb)] . (b) As tension builds 
following contraction, the motors detach and reorganize and 
contractility diminishes. 



reversal symmetry discussed in this paper, and the re- 
sulting requirements for contractility are unchanged. 

A more interesting question is to ask whether contrac- 
tility could arise from the load-dependent detachment 
of myosin motors. Specifically, the detachment rate of 
myosin motors decreases under increasing load, a ten- 
dency known as the "Fenn effect" [31] • To assess its in- 
fluence on contractility, we consider the simple bundle 
of rigid filaments illustrated in Fig. [5](a). This bundle 
comprises a contracting and an expanding region simi- 
lar to those of Fig. [ijb) . When faster motors happen to 
be concentrated in the contracting region, this bundle is 
contractile. Assuming that its ends are fixed, the bundle 
comes under tensile force during contraction. As a conse- 
quence, its expanding region comes under negative load 
and the motors there tend to detach. Conversely, the 
motors in the contracting region experience a positive 
load and tend to hold on to the filaments. On average, 
detached motors thus tend to diffuse from the expanding 
to the contracting region [Fig.^b)]. The reduced num- 
ber of motors in the expanding region implies that its 
resistance to the applied tension decreases, thus increas- 
ing its expansion rate. As the contracting region tends 
to recruit the slower motors from the expanding region, 
its contraction rate decreases. Overall, the Fenn effect 
thus tends to suppress contractile configurations rather 
than amplify them. This suggests that the Fenn effect 
cannot generate contractility in the absence of polarity 
organization or filament nonlinear elastic response. 

The requirements derived in this paper offer insight 
into the contractility of actomyosin bundles devoid of 
sarcomere-like organization, for which no established 
contraction mechanism exists. As they do not satisfy 
condition (2a), we propose that they contract by fulfill- 
ing conditions (1) and (2b) [30| . Consider two antipar- 
allel filaments interacting through several different mo- 
tors with distinct speeds [Fig. [6^ a)]. As motors start to 
move relative to the filaments, stresses build in sections 
of the filament flanked by motors with different speeds. 
When the flanking motor proximal to the barbed end is 
faster than that proximal to the pointed end, compres- 
sion arises. When it is slower, tension arises. F-actin 
responds nonlinearly to these stresses by buckling under 
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Figure 6. Filament buckling as a mechanism for bundle con- 
tractility, (a) The presence of fast (white) and slow (gray) mo- 
tors generically induce compressive (red) and extensile (blue) 
stresses in filaments, (b) Buckling of the compressed filaments 
leads to an overall shortening of the bundle. Buckling of fil- 
ament units of length Li = 500 nm occurs for compressive 
forces of order k B Tt p /L'f ~ 0.16 pN, where £ p ~ 10 /im is the 
F-actin persistence length. Myosin thick filaments typically 
exert forces of several piconewtons [34], which is sufficient 
to induce buckling, but remains smaller than the tension re- 
quired to break F-actin (~ 100 pN [35]). 



compression while resisting extension, which has previ- 



ously been proposed to play a role in actomyosin con- 
traction [32] ■ Following buckling of the compressed fila- 
ment sections, fast motors are free to move quickly while 
the others move slowly. This results in the growth of 
the compressed sections and shrinkage of the extended 
ones, and thus in overall bundle contraction [Fig. |6jb)]. 
Experimental observations suggest that this mechanism 
could be at the origin of contraction in reconstituted ac- 
tomyosin bundles [3U] , and that F-actin buckling can oc- 
cur in cells [33]. These results offer an interesting new 
perspective on mechanisms underlying actomyosin bun- 
dle contraction in vivo. 
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